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(3JT), The time dependent quantum variational principle is emerging as an im- 

portant means of studying quantum dynamics, particularly in early universe 
^ . scenarios. To date all investigations have worked within a Gaussian frame- 

work. Here we present an improved method which is demonstrated to be 
superior to the Gaussian approach and may be naturally extended to the 
field-theoretic case. 



I. INTRODUCTION 

Early Universe scenarios tend to be based on the evolution of scalar fields, either through 
their role as inflaton fields or as topological defect forming fields. A detailed understanding 
of the quantum evolution of these fields is important in order to fully describe their behaviour 
during a phase transition. This has recently been the focus of a great deal of attention. Its 
study relies on analysing the quantum dynamics in real time. Guth and Pi [1] in one of the 
land mark papers in the field of inflation investigate the quantum mechanics of the scalar field 
in the new inflationary universe, looking in detail at the "slow rollover" transition. The idea 
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of a "slow rollover" arises because the transition involves a scalar field which evolves slowly 
down its potential starting from some initial position where it is described by a well defined 
wave function. The phase transition can be thought of one where at very high temperatures 
the potential has a minimum at = 0, which becomes unstable as the temperature decreases, 
with the stable minima moving to a new larger value of = ±o~ say. As the universe cools the 
field remains close to = 0, slowly evolving towards its true vacuum value. The beginning 
of the phase transition is quantum mechanical in nature, yet the late time evolution of 
the scalar field being described by classical equations of motion. This assertion needs to be 
justified, as first addressed in [1]. Analysing an exactly soluble linearised model, both in one- 
dimensional quantum mechanics and in quantum field theory of a single scalar field, it was 
discovered that the large-time behaviour of the field in an unstable upside down harmonic 
oscillator potential is 'accurately described by "classical physics."' In the one-dimensional 
quantum-mechanical model the evolution of the wave function describing the particle in 
the potential is determined by the exactly solvable Schrodinger equation. The solution for 
the wave function is not surprisingly that of a Gaussian. The harmonic oscillator potential 
maintains the form of the initial Gaussian wave function. 

In [2] the time-dependent variational method developed by Jackiw and Kerman [3] is 
used in the investigation of the behaviour of a particle moving in a one dimensional quantum 
mechanical with more realistic potentials (i.e. a double well potential), which have analogues 
both in the inflationary universe scenario and in models for the formation of topological de- 
fects. The important point that needs to be raised, and the motivation for this paper, is 
that in that work, and most subsequent work, the analysis is performed using a Gaussian 
trial wavefunction. The resultant equations of motion obtained through the variation of 
the effective action are the time-dependent Hartree-Fock (HF) equations. In [2] the authors 
argue that by comparison with the exact (numerical) solution the variational HF approxi- 
mation accurately describes the process, and that the late time behaviour of the evolution 
is approximately classical if described in terms of a suitably chosen small dimensionless 
coupling constant. We argue in this paper, that it is straightforward to go beyond this 
Gaussian ansatz by expanding the wavefunction in a complete set of Hermite polynomials. 
In particular we find that just keeping the first and second order terms in the expansion 
leads to a dramatic improvement in the accuracy of the variational approach, and argue 
that such a technique can easily be adapted to the field theory [4] case where presently the 
Hartree- Fock approximation is generally adopted in variational calculations applied to the 
early Universe [5]. 

There are many reasons why it would be advantageous to go beyond the Gaussian ap- 
proximation. If we have a potential which has degenerate minima, then it is impossible 
for a Gaussian wavepacket to accurately describe the evolution of a scalar field during a 
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phase transition. This is particularly relevant for calculations involving the formation of 
topological defects. For example imagine we wish to understand the circumstances under 
which defects can be said to have formed and their distribution at formation [6]. Recently in 
the context of 4 He vortices, there has been considerable attention paid to understanding the 
evolution of the scalar field responsible for their formation just after the quench transition 
[7]. One of the limitations of this interesting calculation is that it can not accurately probe 
the non-linear regions of the potential, hence is only strictly valid just after the quench. In 
order to fully describe the formation process it is important to be able to probe the true 
vacuum of the potential. 

In the context of our simple one dimensional quantum mechanical system we hope to 
demonstrate that by extending the ansatz of the wavefunction to include the Hermite poly- 
nomials, then at little extra cost in complexity we can probe the non-linear region of the 
potential in far greater detail than has previously been possible; an advantage that will be 

carried over to the field-theoretic treatment. 

Another area where the non-linearities of the theory need to be probed, is in the reheat- 
ing calculations associated with inflation. As the inflaton field evolves down its potential, 
eventually it moves out of the "slow roll" regime as it descends into the true minima of 
the potential. The traditional picture is that in this region as the field oscillates about this 
minima then it decays through coherent oscillations and reheats the universe, restoring the 
radiation dominated universe. Recently though, this picture has been questioned [8,9]. In 
order to fully probe this region of the potential, it is important that the ansatz adopted for 
the scalar field is valid in this region. It is our belief that the usual Hartree-Fock approxi- 
mations are not sufficient here and need to be improved. The method we outline in the rest 
of this paper is one possible way of improving the situation. 

Use of the time-dependent variational principle of quantum mechanics is becoming more 
widespread. The method, first given by Dirac, involves the construction of an 'effective 
action' 1 

T = J dt(4>\id t -H\i>), (1) 

where dt = d/dt and H is the Hamiltonian operator. V is then made stationary 5T = 
against variations of the state subject to the constraint (ip\ip) = 1. Approximate 2 



x The connection of which to the usual effective action, or Gibbs free energy is given in [3]. 

2 It is not actually clear in what way the dynamics so obtained are approximate, an issue we will 
address in a later paper. 
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dynamics are obtained by positing a variational ansatz for the wavefunction which is a 

function of a small number of variables. 

Central to the approach is the assumption that the ansatz for the variational wavefunction 

is 'close' to the exact one, i.e. there are sufficient degrees of freedom for the wavefunction 

to accurately track the evolution of the system. 

To date investigations have been restricted to the use of Gaussian ansatzes since these 

are calculationally easy to handle. However we shall argue that the Gaussian approach is of 

limited applicability and that results gained from it have a limited range of reliability. 



II. THE VARIATIONAL PROCEDURE IN GENERALITY 

In order to elucidate our later calculations it is useful to consider this construction in 
generality. Consider the variational effective action 

r = J dtMidt - Hty) = J dt l - «v>IW> - (WIM) - (mm. (2) 

Let us suppose the variational state to be a function of n real parameters Vi 

When the action is made stationary with respect to variation of these parameters we obtain 
the induced equations of motion 



OVj OVi 



v l -l-^\H^) = 0. (4) 



Schematically this expression is of the form 

AijVj - hi = (5) 

implying that in order that we are able to extract the equations of motion for the parameter 
Vj, the matrix A must be non-singular throughout the evolution of the system. 

Obviously thus far we have simply expressed the process one performs in arriving at the 
approximate equations of motion. However, given in this form, the equations of motion are 
simple to extract and the origin of possible singularities is highlighted. 



III. IMPROVED WAVEFUNCTION 

Before we give the expression for the improved wavefunction let us make the definition 
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Unix) := (^—Y H n (ax)e~ a2x2 / 2 a = (2G)"i (6) 

where G(t) is real and H n is the nth Hermite polynomial. We see that the u n are a one 
parameter set of orthonormalized functions 

J u n (x)u m (x)dx = S nm (7) 

independent of the value of G. We may thus use the basis for our variational 

wavef unction: 



oo 

,2 



ifj(x,t) = Ne lUx Y, a nu n (x,t). (8) 

n=0 

Here the a n are time dependent complex numbers, the other real variational parameters being 
IT and G which implicitly appear in the definition of the u n . Normalization is achieved by 
the inclusion of a time dependent N. 

As it stands this represents no simplification. Our plan of action is therefore to truncate 
the expansion at some finite order and work consistently to that order. The zeroth order 
approximation is simply the Gaussian, Hartree-Fock approximation. 

As a test bed for our method we will consider the case of a particle moving in the potential 
considered by Cooper et al. [2], as this illustrates well the shortfalls of the Gaussian approach. 
In a later paper the method will be applied to a wider range of systems. 

The potential we consider is of a double well, 

V(x) = ±(x 2 -ar (9) 

with Gaussian initial conditions Go = ^3/2A a 2 , and a is the symmetry breaking value for 
x. 

To demonstration of the power of the method we will only include the first non-trivial 
term in the expansion. Since the potential and initial conditions are symmetric, the first 
non-trivial term involves u 2 

ifj = Ne iUx2 (u (x, t) + a 2 u 2 (x, t)). (10) 

Even to this order we shall see that the improvement in the results over those of the Gaussian 
approach is impressive. We shall compare the results obtained using the improved equations 
of motion with the exact results obtained via numerical simulation. 
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IV. EQUATIONS OF MOTION 



Working with a polar representation of a 2 (t) = Re ld , the equations of motion one obtains 
from (4) are 

n= 1 _ 2W + ^_JGX_W2cG 

8G 2 12 12 2AR y ' 

it = sAC? (£±^!£±^2±^V2) (13) 

6R 

■ \G 2 UR 3 s/2c+ 2c 2 R 3 - 2c 2 + 1 - 6Ry/2c - 11B 2 ) 1 

9= ulp a (14) 

where s = sin 9 and c = cos#. The increase in complexity of the result over that of the 
Gaussian approach is more than compensated for by the increase in the accuracy of the 
results. 

We notice that the improved equations of motion have within them terms familiar from 
the Gaussian approach. However, the HF equations of motion are not obtainable as a 
simple limit a 2 — > as the improved equations are singular in this limit. The origin of this 
singularity is as was outlined in section (II). 

To assess the use of method we shall focus on the evolution of the quantity 

( .. 2) = gq+^+yggc) {w) 



V. INITIAL CONDITIONS 

It is clear that the equations of motion are singular for the Gaussian initial conditions 
considered; since our aim is a comparison of the results of the improved method with those 
of Cooper et al, we shall adopt the same initial conditions. This apparent problem is easily 
circumvented. 

The rationale is to start the evolution of the system some short time after t — 0. We 
solve the exact Schrodinger equation to first order and may then extract the values of R and 
9, using these as the initial conditions in the variational equations of motion. In actuality 
the subsequent evolution of the system is rather insensitive to the initial conditions. 

Following [2], we take as the initial conditions 

Mx) = (2TrG )^e' x2/4Go (16) 
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where G = 



After a short time St this evolves into 



ip(x, St) = ipo(x) — iStHipQ^x) 



(17) 



which is the state we match onto. This procedure leads us to take 



Rq — 



StV2(6 + Aa 4 G ) 
24 




(18) 



VI. RESULTS AND CONCLUSIONS 



The results of the above calculation are presented in figures VI and VI. Plotted are the 
exact evolution of (x 2 ), found by numerical simulation, against the improved and Gaussian 
results for the two values a = 5,7. It is clear that the improved method furnishes us with 
a result considerably closer to the exact evolution than does the Gaussian. Also we see 
that the improved method samples regions of the potential much closer to the minima than 
does the HF, the so-called spinodal regions, suggesting that more information about the 
potential is being taken into account. The Gaussian wavefunction leads to a turning point 
of (x 2 ) at 2 /3a 2 [2], which provides an indication of where the ansatz breaks down. With 
the improved ansatz we find the turning point occurs typically at (x 2 ) ~ a 2 , demonstrating 
the significant increase in accuracy. Moreover, this approximate solution clearly probes the 
non-linear region of the potential. 

Behind this success is the crucial observation that the improved wavefunctions are capa- 
ble of becoming bimodal in nature, something not open to the Gaussian ansatz. This means 
that we have a method of investigating the field evolution during a defect forming transition 
and a first order transition. These cases are currently being analysed. 

Perhaps the most promising aspect of this work is the possible extension to the field 
theoretic case. We may identify the Hermite polynomials as the eigen-solutions of the 
harmonic oscillator. In field theory the analogous system is the massive free-field. The 
eigensolutions for this system are easily calculable and may then be used as a basis, as 
above. This work is on-going. 
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